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Phase Retrieval 

Jameson Cahill, Peter G. Casazza, John Jasper, and Lindsey M. Woodland 

Abstract. We answer a number of open problems concerning phase 
retrieval and phase retrieval by projections. In particular, one main the¬ 
orem classifies phase retrieval by projections via collections of sequences 
of vectors allowing norm retrieval. Another key result computes the 
minimal number of vectors needed to add to a frame in order for it to 
possess the complement property and hence allow phase retrieval. In 
furthering this idea, in a third main theorem we show that when a col¬ 
lection of subspaces is one subspace short from allowing phase retrieval, 
then any partition of orthonormal bases from these subspaces into two 
sets which fail to span, then each spans a hyperplane. We offer many 
more results in this area as well as provide a large number of examples 
showing the limitations of the theory. 


1. Introduction 

In the setting of frame theory, the concept of phaseless reconstruction 
was introduced in 2006 by Balan, Casazza, and Edidin [I]. At that time, 
they classified phaseless reconstruction in the real case by proving that a 
generic family of (2N-l)-vectors in 1Z N does phaseless reconstruction and no 
set of (2N-2)-vectors can do this. In the complex case, they showed that a 
generic set of (4N-2)-vectors does phaseless reconstruction. Since then there 
has been a multitude of mathematical research devoted to this area. One 
area in particular is the study of phaseless reconstruction by projections 
onto subspaces. An in depth study of phase retrieval by projections was 
done by Cahill, Casazza, Peterson, and Woodland in [6j where they showed 
that phase retrieval by projections can be done in the real case with (2 N — 
l)-projections of arbitrary non-trivial rank and in the complex case with 
(4 N — 2)-projections. We continue the study of phase retrieval by providing 
a variety of new results in both the one-dimensional and higher dimensional 
cases as well as giving a large number of examples showing the limitations 
of the theory. 

The complex case has proven to be a very difficult problem and as such 
has yet to be fully classified. There has been incremental progress towards 
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this end and in particular Heinossaari, Mazzarella and Wolf [9j showed that 
M-vectors doing phaseless reconstruction in C N requires M > AN — 4 — 2a, 
where a is the number of 1/s in the binary expansion of (IV— 1). Bodmann [4] 
later showed that phaseless reconstruction in C N can be done with (41V — 4)- 
vectors. Following this, Conca, Edidin, Hering, and Vinzant [8] proved 
that a generic frame with (41V — 4)-vectors does phaseless reconstruction in 
C N . They also showed that if N = 2 k + 1 then no set of M-vectors with 
M < AN — 4 can do phaseless reconstruction. Bandeira, Cahill, Mixon, 
and Nelson [2] conjectured that for all N, no fewer than (41V — 4)-vectors 
can do phaseless reconstruction. Recently, Vinzant m showed that this 
conjecture does not hold by giving 11 vectors in C 4 which do phase retrieval. 
And recently, Xu m showed that phase retrieval can be done in 7Z 4 with 
six 2-dimensional subspaces. 

2. Frame Theory 

This section first provides a basic review of necessary terms and theorems 
from finite frame theory and then goes on to develop numerous new results 
regarding frames, Riesz bases and complements. For a more in depth study 
of finite frame theory the interested reader is referred to Finite Frames: 
Theory and Applications [7]. To set notation, given N € A f, FLn represents 
a (real or complex) Hilbert space of (finite) dimension N. 

Definition 2.1. A family of vectors {pi}f£i in an N-dimensional Hilbert 
space FLn is a frame if there are constants 0<A<B<ooso that for all 
x € FLm, 

m 

A\\x\\ 2 < \{x,Vi)\ 2 < B\\x\\ 2 , 

1=1 

where A and B are lower and upper frame bounds, respectively. The 
largest A and smallest B satisfying these inequalities are called the optimal 
frame bounds. 

(1) If A = B is possible, then {(pi}f£i is an A-tight frame. If A = 
B = 1 is possible, then {<p>i}££ 1 is a Parseval frame. 

(2) If ||(/?j|| = 1 for all i € [M] then is an unit norm frame. 

(3) {{x, <Pi}}f£i are called the frame coefficients of the vector x € FLn 
with respect to frame 

(4) The analysis operator of the frame is T : FLn —> given by 

T(x) = {(x,^)}^. 

(5) The frame operator S of the frame, for any x € FLn, is given by 

M 

Sx = y ^(x,<pi)^i. 

1=1 

The frame operator is a particularly important object in the study of 
frame theory. 
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Theorem 2.2. Let {ipi}f£i be a frame for LLn with frame bounds A and 
B and frame operator S. S is positive, invertible, and self adjoint. Moreover, 
the optimal frame bounds of are given by A = A m in(S) and B = 

A max(S), the minimum and maximum eigenvalues of S, respectively. 

In particular, the frame operator of a Parseval frame is the identity op¬ 
erator. This fact makes Parseval frames very helpful in applications because 
they possess the property of perfect reconstruction. That is, {-Ai}££i is a 
Parseval frame for LLn if and only if for any x € LLn we have 

M 

x= y 
2=1 

There is a direct method for constructing Parseval frames. For M > N, 
given an M x M unitary matrix, select any N rows from this matrix, then the 
column vectors from these rows form a Parseval frame for LLn ■ Moreover, 
the leftover set of M — N rows also have the property that its M columns 
form a Parseval frame for LLm—N ■ The following well known theorem, known 
as Naimark’s Theorem, utilizes this type of construction and shows that 
this is the only way to obtain Parseval frames. 

Theorem 2.3 (Naimark’s Theorem; page 36 of |7j). Let = {ifi\iL\ be 
a frame for LLn with analysis operator T, let {ei}f£ x be the standard basis 
of and let P : I 2 (M) —>• £2 (M) be the orthogonal projection onto 

range (T). Then the following conditions are equivalent: 

(1) is a Parseval frame for LLn- 

(2) For all i = 1,... , M, we have Pei = Tipi. 

(3) There exist 'if 1 ,... ,ipM £ LLm-n such that {ipi © ifi\fL x is an or¬ 
thonormal basis ofLLM- 

Moreover, if (3) holds, then is a Parseval frame for LLm-n- 

If {V’ili'Ju i s another Parseval frame as in (3), then there exists a unique 
unitary operator L on LLm-n such that Lifi = ipl, for all i = 1,..., M . 

Explicitly, we call the Naimark Complement of <J>. If <J> = 

{<Pi}f£ 1 is a Parseval frame, then the analysis operator T of the frame is an 
isometry. So we can associate ipi with Tipi = Pei, and with a slight abuse 
of notation we have: 

Theorem 2.4 (Naimark’s Theorem). 4> = is a Parseval frame 

for LLn if and only if there is an M-dimensional Hilbert space K-m with an 
orthonormal basis {ej}(£ 1 such that the orthogonal projection P : Km —> LLn 
satisfies Pei = ipi for all i = 1,..., M. Moreover, the Naimark complement 
0 /<h is {(I - P)ei}fi v 
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Note that Naimark complements are only defined for Parseval frames. 
Furthermore, Naimark complements are only defined up to unitary equiva¬ 
lence. That is, if C Hn and C Hm _jv are Naimark com¬ 

plements, and U and V are unitary operators, then {Uipi}f£ 1 and 
are also Naimark complements. 

To clarify terminology, as mentioned in Naimark’s Theorem and as will 
be throughout this paper, an orthogonal projection or simply a projection is 
a self-adjoint projection. 

Frames in the finite dimensional setting are just spanning sets. At times, 
it is useful to look at subsets of a frame which are also spanning sets. 

Definition 2.5. A frame in Hn satisfies the complement 

property if for all subsets S C {1 either span{</?j}j e s = Hn or 

span{^j} ie5 c = U N . 

Definition 2.6. Given a family of vectors in Hn, the 

spark of <b is defined as the cardinality of the smallest linearly dependent 
subset of <b. When spark(<f>) = N + 1, every subset of size N is linearly 
independent and is said to be full spark. 

Remark 2.7. Let <f> := 1 be a frame in Hn such that M > 2N—1. 

(1) If & is full spark then has the complement property. 

(2) If has the complement property and M = 2N — 1 then is full 
spark. 

The notion of spark is the measure of how resilient a frame is against 
erasures, so full spark is a desired property of a frame. In general, it is very 
difficult to check the spark of a frame. In attempts to further characterize 
full spark frames, in [5| the authors classified full spark Parseval frames 
through the use of the frame’s Naimark complement. 

Proposition 2.8. [5] A Parseval frame is full spark if and only if its 
Naimark complement is full spark. 

Often times in applications, a frame is linearly dependent and hence the 
decomposition of a signal with respect to a frame is not unique. However, 
it may be necessary to have a unique decomposition without restricting 
the frame to such properties as orthogonality. A Riesz basis provides this 
uniqueness and does not have as strong of a condition as orthogonality. 

Definition 2.9. A spanning family of vectors {cpi}((L 1 in a Hilbert space 
Hn is called a Riesz basis with lower (respectively, upper j Riesz bounds 
A (respectively, B ), if, for all scalars {aj}-^ 1 , we have 

N N N 

A^T \di\ 2 < II ^a^ill 2 < B^r \ai\ 2 . 
i=l 2=1 2=1 

Recall a result from [5j: 
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Theorem 2.10. Let P be a projection on TLm with orthonormal basis 
{ei}f£ 1 and let I C {1, 2,..., M}. The following are equivalent: 

(1) {Pet} i^i is linearly independent. 

(2) {(/ - P) ei} i£l c spans (/ - P) Ti. 

As a consequence we have, 

Corollary 2.11. Let P be a projection of rank N on TLm with or¬ 
thonormal basis {e^}. The following are equivalent: 

(1) {Pei}ff 1 has the complement property. 

(2) Whenever we partition {(I — P)ei}f£ 1 into two sets, one of them is 
linearly independent. 

The previous results analyze when a collection of projections and their 
orthogonal complements are linearly independent and when they possess the 
complement property. A natural next step is to see when a projection of an 
orthonormal basis is full spark. 

Proposition 2.12. Let P be a projection of rank N on TLm with or¬ 
thonormal basis {ei}f£ 1 . The following are equivalent: 

(1) {Pei}fL l is full spark. 

(2) For every L C {1,2,... M} with |I| = M — N the vectors {(I — 
P) e i} i £ j spans (I — P){TL). 

Proof. This follows immediately from Proposition 12.81 □ 

Instead of only considering an orthonormal basis for our projections, the 
next proposition slightly generalizes this idea by using Riesz bases. 

Proposition 2.13. Let {pi}^L 1 be a Riesz basis with dual basis {p*}^ 
for TLn and let P be an orthogonal projection on TLn of rank M. Let I C 
{1,..., IV}. The following are equivalent: 

(1) {P<Pi} i£i spans PTLn 

(2) {(I-P)p*} i<zic is independent. 

Proof. (1) => (2) (Proof by contrapositive.) Assume {(I — P)p*}i € ic is 
NOT independent. Choose not all zero so that Yliel c bi(I — P)p* = 

0. Then x := Y^iel c hp* = )C,; e /c biPp* € PTLn- If j € I, then 

(x, Pipj) = (Px , tpj) = (x, ipj) = = 0 

i£l c 

since i € I c and j € /. Thus x T span {Ppj}j e j showing that {Ppjfj^j 
does not span PTLn- 

(2) =>■ (1) (Proof by contrapositive.) Assume span {Ppi}^ / PLLjy. 
Then there exists a non-zero x € PTLn with x _L span {Ppifi^j. Also 
x = J2iLi( x > Vi)T*- Now f° r * G I, (x, Ppi) = {Px , pi) = (x, = 0. Hence 

^2(x,pi)p* = x = Px = Y^(x,Pi)Pp*- 

i&I c i£l c 
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Thus, — P)tp* = 0 where (x,<pi) / 0 for at least one i € I c . 

Therefore {(/ — P)(^*}; G /c is NOT independent. □ 

Similarly, the following result classifies full spark projections of Riesz 
bases as follows. 

Proposition 2.14. Let {ipi\f =1 be a Riesz basis on TLn and P an or¬ 
thogonal projection on TLn of rank M. The following are equivalent: 

( 1 ) {Pipi}P =1 is full spark. 

(2) For all I C {1, 2,... , iV} with |/| = M we have: 

[span{(/ - P)ipi}? =1 ] O spanj^Jig/ = { 0 }. 

Proof. {P^Pi}fL 1 is not full spark if and only if there exists an I C 
{1 ? 2,..., TV"! - with |J"| — M and not all zero such that ^ ' jcj arPifii = 0 

if and only if a,i{I ~ P)Vi = a iVi and on ly 

[span{(/ - P)^}^] O spanj^ljg/ ^ { 0 }. 

□ 

A natural question to ask, in light of Proposition 12.141 is if is a 

Riesz basis for TLn and {P^ij^li is full spark on its range for some rank-M 
projection P, then is {(/ — P)'*Pi,}f = \ full spark on its range? The following 
example shows that the answer is no. 

Example 2.15. Let {ei,ei + e 2 } be a Riesz basis for 7Z 2 , where {ei}? =1 
is the standard orthonormal basis for TZ 2 . Let P be the rank-1 projection 
onto e\. Then {Pe\, P(ei + 62 )} = {ei} is full spark on its range. However, 
{(/ — P)ei, (I — P)(e 1 + 62 )} = {0, 62 } is not full spark on its range. 

3. Phase Retrieval 

Phase retrieval has been a long standing problem in mathematics and 
engineering alike. Recently, the mathematical study of phase retrieval by 
subspace components (or projections) has been more deeply developed and 
we further study this area here. In particular, this section expands and 
generalizes a few results from [ 6 ] and |5j as well as develops many new results 
with its main theorem focused on classifying phase retrieval by projections 
via norm retrieval by vectors. 

Definition 3.1. A family of vectors does phase retrieval in 

TL if whenever x,y € P satisfy 

\{x,Vi)\ = \{y,Vi)\, for all i = 1 , 2 ,... ,M, 
then x = cy for some |c| = 1 . 

A family of projections {Pi,}ff \ does phase retrieval in TL if whenever 
x, y £ TL satisfy 

\\Pix\\ = \\Piy\\, for alii = 1,2, 
then x = cy for some |c| = 1 . 
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If{Wi}f£ 1 is a family of subspaces ofH, we say they do phase retrieval 
if the orthogonal projections onto the Wi do phase retrieval. 

One of the main results in [lj is the following. 

Theorem 3.2. A family of vectors $ in R N does phase retrieval if and 
only if <f> has the complement property. In the complex case, phase retrieval 
implies complement property but the converse fails. 

In [6], the authors analyzed phase retrieval by projections and provided 
the following results. 

Theorem 3.3. We have: 

(1) Phase retrieval in IZ N is possible using 2N — 1 
any dimension less than N. 

(2) Phase retrieval in C N is possible using 4 N — 3 
any dimension less than N. 

Also, recall the following result from [5|. 

Proposition 3.4. If is a frame in Hn which allows phase re¬ 

trieval then {Ppi]fLi allows phase retrieval for all orthogonal projections P 
on Hn- 

Generalizing Proposition 13.41 to subspaces leads to the following result. 

Proposition 3.5. If is an orthonormal basis for % 2 N, then for 

all k < N there exists a projection of H. 2 N onto a k-dimensional subspace 
W such that {Pei}^\ does phase retrieval on its range. 

Proof. Without loss of generality assume k = N, choose any Parseval 
frame {<7i }‘i=\ which does phase retrieval for an IV-dimensional Hilbert space 
Kn- Let T be its analysis operator, then T is an isometry and T : Kj\; —>• 
with Tgi = Pei for all i = 1,...,2IV. Thus P is a projection onto the 
range of T, a ^-dimensional space, and {Pe,;}^ = {Tgi}'jil x gives phase 
retrieval. □ 

We can further generalize Proposition 13.51 to Riesz bases as follows. 

Proposition 3.6. If $ = {pi}fL\ is a Riesz basis for %m and N < M 
then there exists a projection P on PLm of rank N so that {Pipi}ff =1 is full 
spark on its range. Hence, if 2N — 1 < M then {P<Pi}f£i yields phase 
retrieval on its range. 

Proof. Let S be the frame operator for <h. Then {S' - is an 
orthonormal basis for %m- Hence, there exists a projection P of rank N 
such that {PS - 2 (^>j}P 1 is full spark. Thus by Proposition [2J3] ( I—P)Hm H 
span{S - 2 <£>j}j e / = {0} for all I C {1,..., M} with |/| = N. 


subspaces, each of 
subspaces, each of 
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Let W = S 2 (/ — P)Pm and suppose x G W n span{(^j}j g /. Then 
x = 5a (I — P)y = Xlie/ for some y G Pm, and bi G iL. So 

S~^x = (/ - P)y = ^ biS~^ipi G (I - n span{5 _ ^j} ie /. 

iei 


Hence 

VL n span{(/?j}j G / = S? (I - P)U M n span{<yj;} ie/ = {0}. 

Now let Q be the projection onto W. Then {(/ — Q)^}^ = {P<p t }ff\ is 
full spark, by Proposition 12.141 When 2N — 1 < M then {(/ — = 

{ Pipi}f£i is full spark and thus has the complement property. Therefore 
{Pipi}^£ 1 yields phase retrieval. □ 


Remark 3.7. The set of vectors {e±,e 2 , e$, ei + e-i, e\ + e^, e^ + 63} have 
the property that it does phase retrieval but no full spark subset does phase 
retrieval. 


Corollary 3.8. If {92, }^ a Riesz basis for Pm with dual Riesz basis 
{<Pi}i£i and 2N — 1 < M < 2N + l, then there exists a projection P of rank 
N so that both of the following hold: 

(1) {Ptpi}f£ 1 does phase retrieval on its range, and 

(2) {(I — P)<p*}fL l does phase retrieval on its range. 

Proof. The proof follows from Proposition 13.61 □ 

In [6], the authors provided necessary and sufficient conditions for when 
subspaces do phase retrieval by relating them to the one dimensional case. 
Their results were proven for the real case. The complex case is more techni¬ 
cal so we will prove this here. To accomplish this, we give a few preliminary 
results. Note that the following results hold in both the real and complex 
case. 


Lemma 3.9. Let be subspaces of Pn allowing phase retrieval. 

For every orthonormal basis ofWi, the set = {Ti,j}i=i)=i allows 

phase retrieval in Pn. 


Proof. Let Pi be the orthogonal projection onto W t for each i = 
{1,..., M}. For every i = 1,2 ,M and for any x,y € Pn such that 
| (a, = \{y, Vi,j)\ for all j = { 1 ,... , Ji}, we have 

Ji Ji 

\\ p ix\\ 2 = Y2\( x ’Vhj)\ 2 = Y2\(y^i,j)\ 2 = \\Piy\\ 2 - 

j =1 3 = 1 

Since allows phase retrieval it follows that x = cy for some c G C 

with |c| = 1. Hence, <f> allows phase retrieval. □ 

Lemma 3.10. Let P be a projection onto an M-dimensional subspace, 
W, of Pn ■ Given x, y G Pn, the following are equivalent: 
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(1) \\Px\\ = \\Py\\. 

(2) There exists an orthonormal basis {pi}f£i for W such that |(x, pf) \ = 
\{y,Vi)\, for all i = 1,2 

Proof. (1) => (2): Consider the vectors Px,Py € W with ||Px|| = 
||Py||. We examine the following three cases. 

Case 1: Assume Px = cPy for some |c| = 1. 

In this case, for any orthonormal basis {pi}f£i for W we have 
\(x,Pi)\ = \(x,Ppf)\ = \ (Px,ipi)\ = \ {cPy,ipi)\ = \c\\(y,Ppi) \ = \(y,Pi)\, 
for all i = 1,2,..., M as desired. 

Hence, for the next two cases, we can assume Px £=■ cPy for any |c| = 1. 
Case 2: Assume ( Px,Py ) = 0. Hence ( Py,Px ) = 0. 

In this case, let 

Px + Py Px — Py 

^ = ||Px + Py\\ ’ and = \\Px — Py\\' 

Note p\ and p 2 are both unit norm. Letting c = l/(||Px + Py||||Pa; — 
Py ||) we have 

(¥> 1 ,^ 2 ) = c(Px+Py,Px-Py) = \\Px\\ 2 — \\Py\\ 2 + (Py,Px) — (Px 1 Py) =0. 

So {p\, P 2 } is an orthonormal set. Also, 

\(x,Px + Py) | = | (Px,Px) + (Px,Py)\ = ||Px|| 2 

= \\ p y\\ 2 = | (Py,Px) + (Py,Py)\ = \{y,Px + Py)\. 

Similarly, \(x,Px-Py)\ = \(y,Px- Py)\. Hence, \(x,pf)\ = \(y,pf}\ for 
i = 1,2. Note that Px,Py £ span {p\, P 2 }- Now, take {pi}f£ 1 to be any 
orthonormal completion of {p\,p 2 \ to an orthonormal basis for W. Then, 
for all i = 3,4,..., M we have 

(x, pf) = (x, Ppf) = (Px, pi) = 0 = (Py, pi) = (y, Ppf) = (y, pf). 

Therefore, \(x,pf)\ = \(y,pf)\ for all i = {1 where {p%}££\ is an 

orthonormal basis for W, as desired. 

Case 3: (Px, Py) / 0. 

In this case, let d = (Px, Py)/\(Px, Py)\ so that |d| = 1, and let 

Px + dPy Px — dPy 

^ = \\Px + dPy\\' an ^ 2 = ||P® - dPy||' 
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Note that and ip 2 are both unit norm. Letting c = 1/1| Px + dPy\\\\Px — 

dPy || we have 


(<Pl,<P2) = 


c(Px + dPy , Px — dPy) 
c(||Px|| 2 - \\dPy\\ 2 + (dPy,Px) - (Px,dPy )) 

c([|Px|| 2 - \d\\\Py\\ 2 ) + d(Py,Px) -d(Px,Py)^j 


C ((l-|<i|)||Fi|| 2 + 

c(o + 


(Px,Py)(Py,Px) _ (Px,Py)(Px,Py) \ 
\(Px, Py)\ \(Px,Py)\ J 

\(Px,Py)\ 2 _ \(Px,Py)\ 2 \ 

\{Px, Py)\ \(Px,Py)\) 

0 . 


Hence, {fi, f 2 } is an orthonormal set. Now the proof follows as in Case 

2 . 


(2) =>- (1): This is immediate by: 

M M 

II Pxf = ^\{x,fi)\ 2 = ^2\(y,f i )\ 2 = \\Py \\ 2 . 

2—1 i=l 

□ 

Combining Lemma [3.91 and Lemma [3.101 we arrive at a characterization 
for when {H 7 j}-^ 1 does phaseless reconstruction in TLn in terms of orthonor¬ 
mal bases. Note, in [0J the authors proved the following result for the real 
case. We will now give the proof for the complex case. 

Theorem 3.11. Let {Wi}f£i be subspaces of TLn- The following are 
equivalent: 

(1) {Wi}^ allows phase retrieval in TLn- 

(2) For every orthonormal basis °fWi, the set i)=i 

allows phase retrieval in TLn- 

Proof. (1) => (2): This is Lemma [3T9l 

(2) (1): Suppose we have x, y € TLn with ||Pjx|| = ||Pjy||, for all 

i = 1,2,... ,M. By Lemma 13.101 we can choose orthonormal bases T = 
so that 

\(x,fij)\ = \ (y,<Pi,j)\, for all i,j. 

By (2), d> does phase retrieval and so x = cy for some |c| = 1. I.e. {Wj}^ 1 
does phase retrieval. □ 

Since orthonormal bases are very restrictive, we would like to relax the 
conditions in Theorem 13.111 to see what properties the vectors within the 
subspaces have when the {Wi}ff 1 are assumed to allow phase retrieval. A 
natural next step would be to look at a Riesz basis as opposed to orthogonal 
vectors. In particular, since unitary operators are the only linear operators 
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which preserve orthogonality, by moving to a Riesz basis we can instead 
look at invertible operators. However, the following example shows that if 
{Wi}fl i allows phase retrieval in 7Z N and i is a Riesz basis for Wi for 

each i 6 then it is not necessarily true that {<Aj}y=;u = i allows 

phase retrieval in 1Z N . 

Example 3.12. Let {ej}f =1 be an orthonormal basis for Tv! 3 . Define the 
subspaces 

W\ = spanjei, e 2 }, W 2 = span{e 2 }, W 3 = span{e 3 }, 
jjr ei + e 2 r e 2 + e 3 f ei + e 3 

Let x € R 3 . Then x = Xa=i a i e i where a t = (x,ei) for i = 1,2,3. We 
have 

af + a I, i = 1 

«2, * = 2 

a|, i = 3 

i(«i + a 2 ) 2 , i = 4 

^ (<T2 T C13) 5 % = 5 

^ («i + 0:3)^ , i = 6 

First, we will show that we can recover ±x from {11-Pppv3?11 2 }f= 1 • 

We can recover the absolute values of the coefficients: 

I«l| = \J\\Pw 1 x\\ 2 - \\Pw 2 x\\ 2 , |CK 2 1 = ||PlV 2 x||, |«3| = ||LV 3 ®||- 

Thus, if two of the coefficients are zero then we have x = ±|aj|ej for 
some i. From now on we will assume at least two of the coefficients are 
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However, if we choose the linearly independent (and not orthonormal) 
basis {ei+e 2 , e 2 } for W\ and the spanning element from the other subspaces, 
we get the set of vectors {e\ + e 2 , e 2 , e 2 , e 3 , e\ + e 2 , e 2 + e 3 , e\ + e 3 } = {e\ + 
C 2 5 e 2 ,e 3 ,e 2 + e 3 ,ei + e 3 }. Notice that if we partition this set as follows 
{^2 + e 3 , e 2 , e 3 }, {e± + e 2 , e\ + e 3 } then neither set spans R 3 . Hence this set 
does not have the complement property and therefore does not allow phase 
retrieval. 

Remark 3.13. In Examvle \3.1A does phase retrieval. To see 

this, we just need to see that this family can retrieve the norm of any vector 
x € 7Z 3 . Let {Qi}f =1 fre the orthogonal projections onto each of {W ^~}® =1 . 
Then, 

W^~ = span{e 3 }, W-f = span{ei,e 3 } = span{ei,e 2 }. 

Given a vector x € IZ 3 , we have 

||(5ix|| 2 = |(z,e 3 }| 2 and \\Q 2 x || 2 = |(a;,ei)| 2 + |(x,e 3 )| 2 

and ||Q 3 x|| 2 = l(^,ei)| 2 + |(x,e 2 )| 2 . 

Hence we know {|(ai, e*)| 2 }^ =1 and so we know ||x|| 2 . 

In light of Example 13.121 we cannot replace “orthonormal bases” with 
“Riesz bases” in Theorem 13.111 (2). The next theorem shows that the key 
property of orthonormal bases in Theorem 13.111 is not just that they span, 
but that they give norm retrieval. 

Definition 3.14. A family of vectors {ipiffL^ does norm retrieval in 
Hn if whenever x, y € Hn satisfy 

\{x,Vi)\ = \(y,Ti)\ for all i = 1,2,... ,M, 

then ||x|| = ||y||. 

A family of projections {Pi}ffi does norm retrieval inH n if whenever 
x,y £ Hn satisfy 

\\Pix\\ = \\Piy\\ for alii = 1 , 2 ,..., M, 

then ||cc|| = ||y||. 

Theorem 3.15. Let {H / j}('£ 1 be subspaces of Hn ■ The following are 
equivalent: 

(1) {H / i})£ 1 allows phase retrieval in Hn- 

(2) For every sequence {Ti,j}j=i C W % which gives norm retrieval in 

Wi, the sequence allows phase retrieval. 

Proof. (1) (2) For each i = 1,..., M let {Ti-.jYjLi be a sequence in 

Wi which gives norm retrieval in Wi. Let x,y gHn such that |(x, (fii.j )| = 

| (y, (pij) | for all j = 1,..., J,, i = 1,..., M. For each i = 1,..., M let Pi be 
the projection onto W t . We have 

\(PiX,ipij) | = \(x,Pmj)\ = \(x,<pij)\ = = \(y,Pi<Pij)\ = \{PiV, Ti,j)\, 
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and since {'•pi.jYjLi gives norm retrieval in W,. we have ||Pjx|| = ||Piy|| for 
all i = 1,..., M. Since {1 / Fj}^£ 1 gives phase retrieval, we have x = cy for 
some c with |c| = 1. 

(2) =>- (1) Since orthonormal bases give norm retrieval, (2) implies that 
each sequence {<Pi,j}j=i i=1 gives phase retrieval, where { l pi,j} J j'=\ is an or- 
thonormal basis for W{. Theorem 13.111 implies that the subspaces 
allow phase retrieval. □ 

Theorem 13.151 provides a clear connection between phase retrieval by 
projections and sequences of vectors giving norm retrieval. In conjunction 
with this theorem, recall the following theorems from [6j and [5|. 

Theorem 3.16. [5J Given i spanning Hn, there exists an invert¬ 

ible operator T on hijy so that the collection of orthogonal projections onto 
the vectors allows norm retrieval. 

Theorem 3.17. |6J If a family of vectors does phase retrieval 

on Hn and T is an invertible operator on Hn, then does phase 

retrieval. 

In light of these results, it seems natural to question if invertible opera¬ 
tors preserve phase retrieval by projections. However, the following example 
illustrates that Theorem 13.171 fails in the higher dimensional case. 

Example 3.18. Let {ej}? =1 be the standard orthonormal basis in 'R?. 
Define the subspaces {H / j}f =1 as in Examvle \3.1S[ Define the linear operator 
T on the basis elements by: 


We have 


Ta 


ei — e2 i = 1 
< e2 i = 2 
,e 3 i = 3 


TWi = W\, TW 2 = W 2 ,TW 3 = W 3 , 

TW A = spanjei}, TW 5 = W 5 , TW 6 = spanjei - e 2 + e 3 }. 

We can choose the following orthonormal bases from each subspace: 


{ei, e 2 } C TWi, {e 2 } C TW 2 , {e 3 } C TW ?J , {ei} C TW A , 
{-/= (e 2 + e 3 )} C TW5, {-^= (ei - e 2 + e 3 )} C TW & . 
We claim that the union of these bases 


{ei, e 2 , e 3 , -^= (e 2 + e 3 ), -^= (e 3 - e 2 + e 3 )} 


does not yield phase retrieval. 

Indeed, if we partition this set into the following two sets: 




1 
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then we see that neither of the two sets spans 1Z 3 . Thus the subspaces 
{TWi}f =1 do not yield phase retrieval. Therefore, we have shown that there 
exists a bounded operator T such that {TWi}® =1 does not yield phase re¬ 
trieval, even though the original subspaces, {Wj}| =1 , did yield phase re¬ 
trieval. 


4. Subspaces which allow phase retrieval 

In this section we develop properties of vectors and subspaces which 
allow phase retrieval. In particular, we start with a frame which fails phase 
retrieval and compute the minimum number of vectors necessary to add to 
the frame so that it will yield phase retrieval. After this computation, we 
then actually construct such sets. The following results show that for a given 
frame which fails the complement property, then for all partitions of the 
frame for which neither set spans 1Z M , there is a minimal number of vectors 
we can add to make this set possess the complement property. Moreover, 
there is an open dense set of vectors which suffice. 

Theorem 4.1. Let be a frame for 1Z N . Fix k £ N L) {0} such 

that 

N — k= min (dim(span{<pj)j g /c)). 

7C{1,...,M} 

spa.n{ip i } ieI ^'R, N 

Then there exists {g>M+ l, • • •, TM+k} C 7Z N such that has the 

complement property. Moreover, if U }y=l h as the complement 

property, then s > k. 

Proof. Let {<pi}f£ l be a frame for 7 Z N . Fix k £ M U {0} such that 

N — k= min (dim(span{<^j}j e /c)). 

/■ {i \i } 

span{^ i } ie /^'R. JV 

Define 


Aq := {I: dim(span{y?j}j e /) < N and dim(span{<^j}j e /c) < N}. 

Since is a finite collection of vectors then there is a finite number 

of sets I Z Aq. For each I C Aq we see that dim(span{<^j}j g /) < N — 1 and 
dim(span{<£>i}, e /c) < N — 1, and hence span{<^j}j e / and span{<^j}j € /c both 
have measure zero. This implies that there is a vector (in fact, an open dense 
set of vectors) pu+i such that tpM+i ^ span{g>j}j e /, tpM+i ^ span{<^j}j € /c, 
dim(span({^j}j e /U <pm+i)) > N-k + 1 and dim(span({^j} ie jc U (Pm+i)) > 

N — k + 1, for any I Z Aq. 

Next, dehne 

A\ := {I: dim(span({^}j 6 /U^M+i)) < N and dim(span({<^} ie /cUv?M+i)) < N}. 
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For each I C A\, dim(span({</?j}j g /U</?M+i)) < IV—1 and dim(span({y>j}j g /cU 
PM+ 1 )) < N — 1 and hence have measure zero. Thus the union of all sub¬ 
spaces defined by / C A\ have measure zero. Hence there exists a vector 
PM+ 2 such that p M + 2 ^ span({</?i} iG / U p M + 1), PM +2 span({</?j} ig jc U 
Pm+i), dim(span({<^}j e /U{^M+j}j=i)) > N-k+2 and dim(span({^i}j g jdJ 
{p M +j}} = i)) > N - k + 2, for any Icii. 

Continue in this manner by defining A^ and adding pM+k until \A^\ = 0. 
Therefore {pi}ff ± k has the complement property. 

For the moreover part let / C {1,... , M} be such that span{<y?j} i&l 7^ 
1Z N and 

N — k = dim(span{y?j}j g /c). 

If {ipj}j= i is a sequence such that {pi}f£i U has the complement 

property, then {pi}i£ic U {ifi} S j=\ must span 1Z N , that is 

N = dim(span({</?j}j g /c U {ipj} S j =1 )) < N - k + s. 

□ 

Theorem 4.2. Let {pt}^ be a frame for 7Z N . Fix k € N U {0} such 
that 

N — k = min (dim(span{(^j} i£l c ))■ 

/e{i,...,M} 

span {ifi lie ipTl N 

Then there exists an open dense set Om+ i C 1Z N such that for all 
Pm+i € Om+i, 

N — k + 1 = min (dim(span{<y9j}j g /c)) . 

/6{1,...,M+1} 

span{(/?i 

Proof. This follows from Theorem O since the union of all subsets 
/ C Aq has measure zero. Hence its complement, call it Om+i-, is full 
measure and thus open and dense. For any vector pm+i € Om+l we have 
PM +1 i- span{<^} ig j, p M +1 £ span {pi} i& i+ dim(span({^} ig / U p M + 1 )) > 
N — k + 1 and dim(span({ 9 Jj}j g /c u pm+ i)) > N — k + 1, for any I <Z A$. □ 

In light of Theorem 14.11 and Theorem 14.21 the next result is surprising. 
It says that if we have families of vectors 4 >i, 4>2 which do phase retrieval 
in two different Hilbert spaces, and we want to do phase retrieval in the 
orthogonal sum of these Hilbert spaces, we will have to add a very large 
number of vectors to 4>i and 4? 2 - 

Theorem 4.3. Let yield phase retrieval on Hnj for j = 

1,2. Let H = "Htvi © 7Ln 2 - If ^ = {Pi)iL\ is a set of vectors in H and 
<f> U 4>i U <f >2 does phase retrieval, then M > N\ + IV 2 — 1. 

Moreover, there exists a 4> = {Pi}^!^ 2-1 so that 4>U4>iU4>2 does phase 
retrieval. 
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Proof. Claim: M > Ni + N 2 — 1 

Proof of claim: If M < N\ + IV 2 — 2, then there exists a partition of 
{1 into {/ 1 ,/ 2 } with |/i| < JV 2 — 1 and |/ 2 | < iVi — 1. Hence 

U doesn’t span for j = 1, 2. 

Proof of the moreover part: Consider the span of all subsets of $1 U $2 
which do NOT span H but do span at least either PLn 1 or 'Hn 2 - There exists 
a vector tpi which is not in the span of any of these subsets. Do the same 
thing with d>iU < I > 2 U{(/?i} and pick (p 2 not in the span of any of these subsets. 
Continue in this way to pick T = 

Claim: <I> U $1 U $2 does phase retrieval. In particular, it has the com¬ 
plement property. 

Proof of claim: Pick a subset J = (Jo, Ji, J 2 ) such that Jo C {1,..., N \+ 
M 2 -l}, J\ C {1,.. •, Mi}, and J 2 C (1,... ,N 2 }. Either {ipj} lG j or {^} ie jc 
spans Hnj for j = 1,2. 

Case 1: If spans T-Lnj for j = 1,2, then we are done. 

Case 2: If spans Pn :i for j = 1, 2, then we are done. 

Case 3: Without loss of generality assume 

U U {$• }; e j 2 spans H Nl 

and {<fo}i G j 0 c U {T’Kej- U {d>f} ieJ c spans Un 2 . 

We have either | Jq\ > iV 2 or | Jf\ > Ni, otherwise 

\J 0 \ + \JS\<N 1 -1 + N 2 -1 = N 1 + N 2 -2, 
a contradiction. 

Without loss of generality, assume that | Jo| > 1V 2 . 

We need to show that {<hj}i e j 0 i S j 2 spans H. 

Let Jo = {k\ < L ’2 < • • • < kjy 2 }. Then 

dim(span({$}} ieJl ,{$ fc J^)) = 1 + dim(span({$}} ieJl , 

= 2 + dim(span({$}} ieJl , {$ fci }^ 2 “ 2 )) 


= 1V 2 + dim(span PL m 1 ) 

= dim('H). 

Therefore $U $1 U $ 2 does phase retrieval. □ 

5. Properties of subspaces which fail phase retrieval 

In the one-dimensional real case, the complement property completely 
classifies phase retrieval. However, the complement property is not a suffi¬ 
cient condition for a collection of subspaces to allow phase retrieval. When 
the complement property fails, we will see that the corresponding partition 
yields two hyper planes. 
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Theorem 5.1. Let {Wi}^L x be subspaces of TLn which yield phase re¬ 
trieval but {Wi}i<zi fails phase retrieval for any I C {1,... , M} with |/| = 
M — 1. If {fi,jYj' =1 is an orthonormal basis for W{ for all i G I, and I\ and 
I 2 is a partition of {(i, j)}^Li so that {fi,j}(i,j)eii f or Z = 1 , 2 do not span 
TLn, then 

dim (span{/ij} (ij)e/l ) = dim (span{/ij} (i>J -) e72 ) = N - 1. 

Proof. Without loss of generality, assume I = — 1}. Let 

Ki = span {fi,j}(ij)ai for t = 

Case 1: K\ ^ K 2 . 

By way of contradiction assume that dim K 1 7 ^ dim K 2 ■ Hence for some 
n > 1, codim(span K\ ) = n and codim(span K 2 ) > n + 1. Then 

dim (Wm H K\) > dim (Wm) ~ n. 

Set W' M := {the orthogonal complement of Wm H K\ inside of Wm}- So 
dim(lT{y-) < n. Next, pick an orthonormal basis Q = {. 9 i}™i for IT^ and 
an orthonormal basis £ = {hi} p i=1 for Wm H K\. Note that m < n. Now, 
QU C is an orthonormal basis for Wm- Add {.9i}™h l 1 f° an( ^ 

9m. to {fij}(i,j)ei 2 - Then neither sets span the space and so {TPj}^ does 
not do phase retrieval - a contradiction. 

Case 2: K\ = K 2 < N — 1. 

By the argument of Case 1, neither of our sets span (because K\ < N— 2) 
and so {H / j }^ 1 does not do phase retrieval - a contradiction. □ 

As a partial converse to Theorem 15. 11 we have the following proposition 
for the one-dimensional case. 

Proposition 5.2. Assume are vectors inH n with the following 

property: 

Whenever we partition {fi}ff 1 into two non-spanning subsets { fi}iGh, 
{fi}iei 2 then each of these sets of vectors spans a hyperplane. 

Then there is an open dense set of vectors fo € TLn so that {fi}f£ 0 does 
phase retrieval. 

Proof. Choose any vector from the open dense set of vectors fo which 
is not in span({/j}j G j) whenever this family does not span the space. If we 
partition {/i}^L 0 i n t° {fi}ieTj, j = 1, 2 , then one of these sets must span the 
space. I.e. If neither set spans the space, then by removing the vector fo 
each family of vectors spans a hyperplane and one of them must contain fo 
and hence spans the space. □ 
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The following corollary is a rephrasing of Theorem 15.11 in such a way 
which may be more useful when proving a collection of subspaces fail phase 
retrieval. 

Corollary 5.3. Suppose {Wi}f £ 1 are subspaces in 7 ~Ln■ Let { ! -Pi 3 }j'L 1 
be an orthonormal basis for Wi for each i = 1 and suppose there 

exists a partition of into two non-spanning sets Ai, A2 • If 

dim(spanFi) < M — 2 then for all subspaces Wm+i of Hn, {Wj}^^ 1 fails 
phase retrieval. 

We would like to know if the converse of Corollary 15.31 is true. This is 
explicitly stated in the following problem. 

Problem 5.4. Let {Vp}^ be subspaces in Fn- If for all subspaces 
Wm+i of Fn , {W^L'f 1 fails phase retrieval then does there exist an or¬ 
thonormal basis {tyijYj-x of Wi and a partition F \, F 2 of {'+ 13 } f=\ );=1 such 
that dim(spanf ? i) < M — 2 and dim(spani ? 2) < M — 1? 

In light of Theorem 15.11 we now further analyze the two hyperplanes 
spanned by a partition of orthonormal bases for our subspaces which fail 
phase retrieval and find properties of the vectors within these hyperplanes. 

Proposition 5.5. Let {Wj}^//" 1 in Fn yield phase retrieval. Assume 
{ptj}j=i an orthonormal basis for Wi, for i = 1 ,..., M. Assume there ex¬ 
ists a partition I\,l 2 of {pij}^Ai\ =l so that dim(span(/ s )) = N — 1, s = 1, 2. 
Then there exists an orthonormal basis {ifij}^L 1 for Wi, for i = 1 ,M 
and a partition Ji, J2 °.f satisfying: 

(1) dim(span( J s )) = N — 1 for s = 1,2, and 

(2) Ji contains at most one vector from i 1 for each i = 1,..., M. 

Proof. Let K s = span(I s ) for s = 1,2. Without loss of generality we 
can assume that there does not exist any tpij € I\ in K 2 or else we can 
move these vectors without affecting K\, K 2 ■ (Note that for any vector cpij 
in 11 which lies in K 2 , we can move this vector to I 2 without affecting K 2 . 
Also, this will not decrease dim(Ad) because if it did then we would be left 
with two non-spanning sets with dimensions N — 2 and N — 1 respectively, 
contradicting the fact that they both must span hyper planes.). So K s is a 
hyperplane for s = 1,2. 

Case 1: If I\ contains at most one vector from for each i = 

{1 ,... ,M} then we are done. 

Case 2: Fix i. Assume I± contains at least two vectors from {<Pij}jf i (so 
dim(span(ITj)) > 2). Let W[ = spanj^jl^jj E I\ for some j E {1,... , Li}}. 
Now codim(span M/ .(A'2 n IP/)) = 1. Choose an orthonormal basis {il>ij}j Z = i 
for K 2 n W[ and choose i/yt ?; +i -L (A 2 n W[) and i/yt ?; +i € IP/. Throw away 
all {<fij}jL 1 which lie in IP/. Put 'fi^+i into I\ and put the rest of {il>ij}j Z =1 
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into I 2 • Repeat this for each i £ {1,..., M}. Define these new sets to be 
Ji, J 2 respectively, and we are done. □ 

Corollary 5.6. For each n £ N let be a sequence in Hn- If 

dim(span{y>J*}jL 1 ) = N — 1 

for all n € N and ipf —> <pi for each i G {1, ..., k}, then 

dim(span{<£j}k =1 ) < N — 1. 
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